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We study the behavior of the gluon and ghost dressing functions in SU(3) Landau gauge at low 
momenta available on lattice sizes 12 4 - 32 4 at p = 5.8, 6.0 and 6.2. We demonstrate the ghost 
dressing function to be systematically dependent on the choice of Gribov copies, while the influence 
on the gluon dressing function is not resolvable. The running coupling given in terms of these 
functions is found to be decreasing for momenta q < 0.6 GeV. We study also effects of the finite 
volume and of the lattice discretization. 
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I. INTRODUCTION 

Studying the relation of non-perturbative features of 
QCD, such as confinement and dynamical chiral symme- 
try breaking, to the properties of propagators, there are 
two popular approaches at present: lattice gauge the- 
ory and Dyson-Schwinger equations (DSE). The latter 
approach allows to address directly the low-momentum 
region for a coupled system of quark, gluon and ghost 
propagators which is of interest for hadron physics 0. 
In particular their infrared behavior could be related to 
the mechanism of dynamical symmetry breaking and to 
confined gluons 0, Q ■ 

In fact, the DSE approach has revealed that in the in- 
frared momentum region a diverging ghost propagator G 
is intimately connected with a suppressed gluon propa- 
gator D^ v . In Landau gauge, they can be written as 



G(q) 



<7 M q v \ Z D {q 2 ) 



Z G {q 2 ) 



(1) 
(2) 



Here Zjj(q 2 ) and Zc(q 2 ) denote the dressing functions of 
the corresponding propagators. They describe the devi- 
ation from the momentum dependence of the free prop- 
agators. Based on the Dyson-Schwinger approach and 
under mild assumptions these functions are predicted to 
behave in the limit q 2 — * as follows pj: 



Z D (q 2 ) <x (q 2 



Z G (q 2 ) oc (q 2 )' 



(3) 



with exponents satisfying kd = 2kg- In Landau gauge 
kq ~ 0.595 Thus the ghost propagator diverges 

stronger than 1/q 2 and the gluon propagator is infrared 
suppressed. This is in agreement with the Zwanziger- 
Gribov horizon condition 0, 0, as we ^ as with the 
Kugo-Ojima confinement criterion Zwanziger |]| has 
suggested that in the continuum this behavior of the 
propagators in Landau gauge results from the restric- 



tion of the gauge fields to the Gribov region f2, where 
the Faddeev-Popov operator is non-negative. 

Using further the ghost-ghost-gluon vertex, the gluon 
and ghost dressing functions can be used to determine 
a renormalization group invariant running coupling in a 
momentum subtraction scheme as 0, Hoi Hll 



a s (q 2 ) 



- ~ 9 lz 2 G ( q 2 )Z D {q 2 ) 



(4) 



which then enters the quark DSE 0, Q|. This defini- 
tion relies on the fact that the ghost-ghost-gluon ver- 
tex renormalization function Zi(q 2 ) is constant, which 
is true at least to all orders in perturbation theory [12j . 
Indeed, a recent numerical investigation of Z\ for the 
SU(2) case shows that also nonperturbatively Z\ is fi- 
nite and constant [lflj . Applying the behavior given in 
Eq. the running coupling approaches a finite value 
a s (0) Ri 8.915/iV for SU(N) at zero momentum in the 
DSE approach 0. 

Nevertheless, numerical investigations of those features 
in lattice simulations are still necessary to check to what 
extent the truncation of the coupled set of DSEs influ- 
ences the final result. There are several studies in Landau 
gauge which confirm the antici pate d behavior at least 
for the SU(2) case El 13 El 1H Also lattice stud- 
ies (0,0,0, |20L l2l| and references therein) for the 
SU(3) case indicate the correctness of the proposed in- 
frared behavior. However, as recent DSE investigations 
show [2^, |3 the infrared behavior of the gluon and 
ghost dressing functions and of the running coupling is 
changed on a torus. In particular, the running coupling 
decreases at low momenta. 

This paper presents a lattice study of the gluon and 
ghost dressing functions and of the running coupling at 
low momenta in SU(3) Landau gauge. We also focus, 
more carefully than usually, on the problem of the Gribov 
ambiguity in lattice simulations. In continuum, a gauge 
orbit has more than one intersection (Gribov copies \J\) 
with the transversality plane (where d^A^ — holds for 
the gauge potential A p ) within the Gribov region f2. Ex- 
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pectation values taken over this region are argued to be 
equal to those over the fundamental modular region A 
which includes only the absolute maximum of the gauge 
functional [5j- 

On a finite lattice, however, this equality cannot be ex- 
pected Q • In the literature it is widely taken for granted 
that the gluon propagator does not depend on the choice 
of Gribov copy, while an impact on the ST/ (2) ghost prop- 
agator has been observed I2M l26l |27| . However, in a 
more recent investigation |28j an influence of the Gri- 
bov copies ambiguity on the St/(3) gluon propagator has 
been demonstrated, too. Here we assess the importance 
of the Gribov ambiguity on a finite lattice for the SU (3) 
ghost propagators. 

This paper is structured as follows: In Sec. [H] we shall 
define all quantities which are investigated in this study. 
Then, after specifying the lattice setup used, the depen- 
dence of the gluon and ghost propagator on the choice of 
Gribov copies and lattice discretization as well as finite- 
volume effects are discussed in Sec. IIIII We also discuss 
the problem of exceptional gauge copies in Sec. I1VI In 
Sec.^the behavior of the running coupling is presented. 
In Appendix ^ we show how the inversion of the F-P 
operator can be accelerated by pre-conditioning with a 
Laplacian operator. 

II. DEFINITIONS 

To study the ghost and gluon propagators using lattice 
simulations one has to fix the gauge for each thermalized 
SU(3) gauge field configuration U — {U x ^}. We adopted 
the Landau gauge condition which can be implemented 
by searching for a gauge transformation 

Ux,ii ► g U x ,^ = g x U x ^ g^fr 
which maximizes the Landau gauge functional 
1 4 

while keeping the Monte Carlo configuration U fixed. 
Here g x are elements of SC/(3). 

The functional Fjj [g] has many different local maxima 
which can be reached by inequivalent gauge transforma- 
tions g, the number of which increases with the lattice 
size. As the inverse coupling constant f3 is decreased, in- 
creasingly more of those maxima become accessible by 
an iterative gauge fixing process starting from a given 
(random) gauge transformation g x . The different gauge 
copies corresponding to those maxima are called Gribov 
copies, due to their resemblance to the Gribov ambiguity 
in the continuum 0. All Gribov copies { 9 U} belong to 
the same gauge orbit created by the Monte Carlo config- 
uration U and satisfy the differential Landau gauge con- 
dition (lattice transversality condition) (d ll 9 A IJ )(x) = 
where 

(d^)(x) = <A^(x + A/2) - "A M (x - A/2)- (6) 



Here 34 M (a; + A/2) is the non-Abelian (hermitian) lattice 
gauge potential which may be defined at the midpoint of 
a link 

«,(x + A/2) = JL.{'U^->>UU 

In this way it is accurate to 0(a 2 ). The bare gauge 
coupling go is related to the inverse lattice coupling via 
f3 — 6/pq in the case of SU(3). In the following, we will 
drop the label g for convenience, i.e. we consider U to 
be already put into the Landau gauge such that g = 1 
maximizes the functional in Eq. (jSJ) relative to the neigh- 
borhood of the identity. 

The gluon propagator D^Kq 2 ) is the Fourier trans- 
form of the gluon two-point function, i.e. the expectation 
value 

DfM = (A;(k)A b u (-k)) = 5 ab D, v {q) (8) 

which is required to be color-diagonal. Here A®(k) is 
the Fourier transform of A®(x + A/2) and q denotes the 
momentum 



= 




which corresponds to a integer-valued lattice momen- 
tum k. Since G (—La/2, Ln/2] the lattice equiva- 
lent of Q 2 (k) can be realized by different k. According to 
Ref. [23, however, a subset of lattice momenta k has been 
considered only for the final analysis of the gluon propa- 
gator, although the FFT algorithm provides us with all 
lattice momenta. Details are given below. 

Assuming reality and rotational invariance we envisage 
for the (continuum) gluon propagator the general tensor 
structure: 

D^(q) = \6„ u - —5- J D(q ) + — (10) 

with D(q 2 ) and F(q 2 ) being scalar functions. On the lat- 
tice these functions are extracted by projection and are 
expected to scatter, rather than being smooth functions 
of q 2 . Using the Landau gauge condition the longitudi- 
nal form factor F(p 2 ) vanishes. Recalling the mentioned 
Gribov ambiguity of the chosen gauge copy there is no 
a priori reason to assume the estimator of D(q 2 ) is not 
influenced by the choice. 

The ghost propagator is derived from the Faddeev- 
Popov (F-P) operator, the Hessian of the gauge func- 
tional given in Eq. (|SJ). We expect that the properties 
of the F-P operator differ for the different maxima of 
the functional (Gribov copies). This should have conse- 
quences for the ghost propagator as is shown below. 

After some algebra the F-P operator can be written in 
terms of the (gauge-fixed) link variables U Xifl as 

^xy = E ^ X >V ~ ^ X +P->V ~ ^a:,p $E~A,V (H) 

A* 
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No. 


13 


lattice 


# conf 


# copies 


selected k for G(k) 






1 


5.8 


24 4 


40 


30 


([1,0], 0,0) 






2 


6.2 


12 4 


150 


20 


([1,0], 0,0) 






3 


6.2 


16 4 


100 


30 


([1,0], 0,0) 






4 


6.2 


24 4 


35 


30 


([1,0], 0,0) 






5 


5.8 


16 4 


40 


30 


"j^f/c^ k ^ k . k^ j k — 1 . . 


•6}, 


(2,1,1,1) 


6 


6.0 


16 4 


40 


30 


k . k . k^ j k — 1 . . 


•6}, 


(2,1,1,1) 


7 


6.2 


16 4 


40 


30 


"j^/Cj k . k . fe) , /ti — 1 . . 


.6}, 


(2,1,1,1) 


8 


5.8 


24 4 


25 


40 


, & ^ k^ i k — 1 . . 


•6}, 


(2,1,1,1), ([1,0,0,0]), (1,1,0,0), (1,1,1,0) 


9 


6.0 


24 4 


30 


40 


"j^/cj k j /c , k^ i k — 1 . . 


•6}, 


(2,1,1,1), ([1,0,0,0]), (1,1,0,0), (1,1,1,0) 


10 


6.2 


24 4 


30 


40 


"j^/cj /c , /i)^ k^ ^ k — 1 . . 


•6}, 


(2,1,1,1), ([1,0,0,0]), (1,1,0,0), (1,1,1,0) 


11 


5.8 


32 4 


14 


10 


(1,0,0,0), (1,1,0,0), (1,1,1,0), (1,1,1,1), (2,1,1,1) 



TABLE I: Statistics of the data used in our final analysis. The 6th column lists all tuples of momentum k the ghost propagator 
has been calculated for. If entries are given in squared bracket all of their permutations are meant. 



with 



A ab 

r>ab 

Cab 



SReTr [{T a ,T b }(U Xtfl 
2-mtTi[T b T a U x ^] , 
2-mtTr[T a T b U x -^] 



Here T a and T b are the (hermitian) generators of the 
su(3) Lie algebra satisfying Tr[T a T b ] = 5 ab /2. 

The ghost propagator is calculated as the following en- 
semble average 



-2-ni k-(x-y) 



[AT 



-liab 
\xy 



(12) 



It is diag onal in color space: G ab (q) = S ab G(q). Following 
Ref. [23, |3(J we have used the conjugate gradient (CG) 
algorithm to invert M on a plane wave ij) c with color and 
position components ip2( x ) = 5 ac exp(2Tri k-x). In fact, 
we applied the pre-conditioned CG algorithm (PCG) to 
solve M^ b c/) b (y) = i/j® (x) . As pre-conditioning matrix we 
used the inverse Laplacian operator A -1 with diagonal 
color substructure. This significantly reduces the amount 
of computing time as it is discussed in more detail in 
Appendix ^ 

After solving M<fi = ip c the resulting vector <j> is pro- 
jected back on ip c such that the average G cc (q) over the 
color index c can be taken explicitly. Since the F-P op- 
erator M is singular if acting on constant modes, only 
k 7^ (0,0,0,0) is permitted. Due to high computational 
requirements to invert the F-P operator for each k, sep- 
arately, the estimator on a single, gauge-fixed configura- 
tion is evaluated only for a preselected set of momenta 
k. In Table [I] a detailed list is given. 



III. RESULTS FOR THE GHOST AND GLUON 
PROPAGATORS 

A. Lattice samples 

For the purpose of this study we have analyzed pure 
SU(3) gauge configurations which have been thermal- 
ized with the standard Wilson action at three values of 
the inverse coupling constant j3 = 5.8, 6.0 and 6.2. For 
thermalization an update cycle of one heatbath and four 
micro-canonical over-relaxation steps was used. As lat- 
tice sizes we used 16 4 , 24 4 and 32 4 . For tests, exposing 
the inherent problems of the Gribov problem under the 
aspect of volume dependence, also smaller lattices (8 4 
and 12 4 ) have been considered at lower cost. 

To each thermalized configuration U a random set of 
N cp local gauge transformation {g x } were assigned. Each 
of those served as starting point for a gauge fixing pro- 
cedure for which we used standard over-relaxation with 
over-relaxation parameter tuned to a; = 1.63. Keeping all 
U x ^ fixed this iterative procedure creates a sequence of 
local gauge transformations g x at sites x with increasing 
values of the gauge functional (Eq. (J5J). Thus, the final 
Landau gauge is iteratively approximated until the stop- 
ping criterion in terms of the transversality (see Eq. ©) 



-14 



(13) 



is fulfilled. Consequently, each random start results in its 
own local maximum of the gauge functional. Certain ex- 
trema of the functional are found multiple times. In fact, 
this happened frequently on the small lattices, 8 4 and 
12 4 , but rather seldom on larger lattices. Note that we 
used the maximum in relation (|13|) which is very conser- 
vative. However, the precision of transversality dictates 
how symmetric the F-P operator M can be considered. 
This is crucial for its inversion and thus dictates the final 
precision of the ghost propagator. 



4 



8.0 - 
7.6 -. 
7.2 - 



12 4 p = 6.2 



k=([l,0],(>,0) I • — I 



-~i - 



14 - 
Q 13 



k=([l,0,0,0]) I — — 



12 



& o 



1(1 

N 



in 



2(1 



13.6 
13.2 
12.8 

2S 



Li 



16 4 p = 6.2 



k=([l,0],0,O) I — • — I 



k=([l,0,0,0]) I — — 

"I i I I I S 



10 15 20 25 30 



24 4 p = 6.2 

' 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
k=([l,0],0,0) h--» 

* + T- 



™ 44^ 

\ 43 ~i 

T\ 42 - r 

: 41 - 



24 4 p = 5.8 



k=([l,0],O,O) I — *-H : 



5 I 



k=([l,0,0,C 



26 r 
25 - 



k=([l,0,0,0]) 



^ i i r xt I 24 5 1 i pit f 

: E 23 - 



10 15 20 25 30 

ff 



2 ^ 
1 =- 
- 




10 15 20 25 30 

N 



FIG. 1: The upper panels show the ghost propagator G(k) as average over two realizations k = (1,0,0,0) and k = (0, 1,0,0) 
of the smallest lattice momentum, measured always on the best gauge copy among N cp copies. In the middle panels the same 
dependence is shown for the gluon D(k) propagator, however, as average over all four permutations of k = (1,0,0,0). The 
lower panels show the relative difference SF = 1 — F cbc / F hc of the corresponding current best functional values 7 ?cbc to the 
value F bc of the overall best copy. 



To study the dependence on Gribov copies of the prop- 
agators, in the course of N cp repetitions for each U, the 
gauge copy with the largest functional value is stored un- 
der the name best copy (be). The first gauge copy is also 
stored, labeled as first copy (f c). However, it is as good 
as any other arbitrarily selected gauge copy. 

The more gauge copies one gets to inspect, the bigger 
the likeliness that the copy labeled as be actually repre- 
sents the absolute maximum of the functional in Eq. JSJ • 
With increasing number iV cp of copies the expectation 
value of gauge variant quantities, evaluated on be repre- 
sentatives, is converging more or less rapidly as we will 
discuss next. 



B. How severe is the lattice Gribov problem for 
the propagators? 

First we present results of a combined study of the 
gluon and ghost propagators on the same sets of f c and 
be representatives of our thermalized gauge field config- 
urations. This allows us to assess the importance of the 
Gribov copy problem. 

Numerically, it turns out that the dependence of the 
ghost propagator on the choice of the best copy is most 
severe for the smallest momentum. In addition, this de- 
pends on the lattice size and (3. Therefore we studied 
first the dependence of the ghost and gluon propagators 
at lowest momentum on the (same) best copies as func- 
tion of the number of gauge copies Ar cp under inspection. 
This was done at (3 = 6.2 where we used 12 4 , 16 4 and 
24 4 lattices. The number of thermalized configurations 
used for these three lattice sizes are given in Table H] To 



check the dependence on j3 also a simulation at j3 = 5.8 
on a 24 4 lattice was performed. 

The results of this investigation are show in Fig. ^ 
While there the ghost propagator is shown as an 
average over the two realizations A; = (1,0,0,0) and 
k= (0,1,0,0) of the smallest lattice momentum q(k), 
the gluon propagator has been averaged over all four 
non-equivalent realizations. Note that D(k) = D(—k). 
It is clearly visible that the expectation value of the 
gluon propagator does not change within errors as A^ cp 
increases, independent of the lattice size and p. Con- 
trarily, the ghost propagator at (3 = 5.8 on a 24 4 lattice 
saturates (on average) if calculated on the best among 
A^cp = 15 gauge copies. At (3 = 6.2 the number of gauge 
fixings attempts reduces to 5 < iV cp < 10 on a 16 4 and 
24 4 lattice. On the 12 4 lattice a small impact of Gribov 
copies is visible, namely 1 < iV cp < 5. The lower panels 
of Fig^show the relative difference 5F = 1— F chc /F hc f 
the corresponding (current best) functional values F cbc 
to the value F bc of the overall best copy after N cp = 20, 
respectively iV cp = 30, attempts. This may serve as an 
indicator how large N cp has to be on average for the cho- 
sen algorithm to have found a maximum of F close to 
the global one. 

In order to study further the low-momentum depen- 
dence of the gluon and ghost dressing functions, Zd and 
Zg, as given by Eq. and Eq. Q we have performed 
similar simulations using lattice sizes 12 4 , 16 4 , 24 4 and 
32 4 at = 5.8, 6.0 and 6.2. Following Ref. ^ these val- 
ues of (3 correspond to a _1 =1.446 GeV, 2.118 GeV and 
2.914 GeV using the Sommer scale tq = 0.5 fm. These 
values of the lattice spacing a associated to (3 turn out 
to be more appropriate as those formerly used by us and 
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FIG. 2: The upper panels show the dressing functions of the ghost Zg and gluon propagator Zd measured on best gauge copies 
as functions of the momentum q 2 (scaled to physical units at f3 = 5.8, 6.0 and 6.2) using various lattice sizes. The lower panels 
show the ratio (Z c )/(Z bc ) determined on first (fc) and best (be) gauge copies using jackknife. 



others (see 

We have fixed a conservative number of N cp = 30 
gauge copies per thermalized configuration on a 16 4 lat- 
tice and N cp — 40 on a 24 4 lattice. Both the gluon and 
the ghost propagator, respectively their dressing func- 
tions, have been measured on the same set of f c and 
be copies. Due to the large amount of computing time 
necessary for the 32 4 lattice we could afford to measure 
the ghost propagator for the first and best among only 
N cp — 10 copies, which is certainly not enough. 

The data for the gluon propagator D(k) have been de- 
termined for all momenta at once. However, we used 
only a subset of momenta for the final analysis. In fact, 
inspired by Ref. |2^, only data D{k) with k lying in a 
cylinder with radius of one momentum unit along one 
of the diagonals n = 1/2(±1, ±1, ±1, ±1) have been se- 
lected. Since we are using a symmetric lattice structure 
only data with k satisfying ^2 k 2 — k^n^) 2 < 1 are 
surviving this cylindrical cut. In agreement with |29| this 
recipe has drastically reduced lattice artifacts, in partic- 
ular for large momenta. Additionally, we try to keep 
finite volume artifacts at lower momenta under control 
by removing all data D(k) with one or more vanishing 
momentum components fc M p9j . However, this we have 
done only for data on a 12 4 and 16 4 lattice. In Sec. 1111 CI 
we shall discuss in more detail finite volume effects at 
various momenta. 

In view of this we have chosen appropriate sets of mo- 
menta for the ghost propagator, as listed in Table [I] in 
detail. 

The final results of the dressing functions Zd and Zq 
measured on be copies are shown in the upper panels of 
Fig. [5] All momenta q(k) have been mapped to physical 
momenta using the lattice spacings a given above. As ex- 



pected the ghost dressing function diverges with decreas- 
ing momenta, while the gluon dressing function decreases 
after passing a turnover at about q 2 = 0.8 GeV 2 . How- 
ever, for the purpose of the expected infrared behavior 
given in Eq. (J3J the data for momenta q 2 < 0.25 GeV 2 
are not sufficiently abundant to extract a critical expo- 
nent kg > 0.5 as expected from the Dyson- Schwinger 
approach. In particular, the fit parameters are not sta- 
ble under a change of the upper momentum cutoff. The 
best fits give kq k 0.25. 

In the lower panels of Fig. [21 we present the ratio of the 
dressing functions, (Z ic )/(Z bc ), calculated using jack- 
knife from first and best gauge copies as a function of 
the momentum. There the data from simulations on a 
32 4 lattice have been excluded, since only A cp = 10 gauge 
copies have been inspected which would result in an un- 
derestimate of the ratio (Z lc )/(Z bc ). As is clear from 
these panels we do not observe a systematic dependence 
on the choice of Gribov copies for the gluon propagator. 
In contrast the ghost propagator is systematically overes- 
timated for f c (arbitrary) gauge copies. This effect holds 
up to momenta of about 2 GeV 2 . 

Comparing also the ratios for the ghost propagator at 
q < 1 GeV, the rise at (3 = 6.0 is obviously larger than 
that at f3 = 5.8. In both cases the data are from simula- 
tions on a 24 4 lattice. Thus, it seems that by increasing 
the physical volume (lower (3) the effect of the Gribov 
ambiguity gets smaller if the same physical momentum 
is considered. This cannot be due to a too small number 
A cp of inspected gauge copies since, judged from Fig. ^ 
N cp = 40 seems to be on the safe side. 

We conclude: the ghost propagator is systematically 
dependent on the choice of Gribov copies, while the im- 
pact on the gluon propagator is not resolvable within our 
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statistics. However, there are indications that the depen- 
dence on Gribov copies decreases with increasing physical 
volume. This is also in agreement with the data listed in 
the two lattice studies (2^, H(j of the SU (2) ghost prop- 
agator G, while it is not explicitly stated there. In fact, 
in Ref. the ratio G lc /G bc at (3 = 2.2 on a 8 4 lattice 
is larger than that on a 16 4 lattice at the same physical 
momentum. 



Systematic effects of lattice spacings and 
volumes 



We remind that in Fig. [21 we have dropped all data 
related to a 16 lattice with one or more vanishing mo- 
mentum components k^. According to |29| this keeps 
finite volume effects for the gluon propagator under con- 
trol. It is quite natural to analyze here the different sys- 
tematic effects on the gluon and ghost propagators of 
changing either the lattice spacing a or the physical vol- 
ume V. However, due to the preselected set of momenta 
for the ghost propagator and the values chosen for f3, we 
can study this here only in a limited way and in a re- 
gion of intermediate momenta. For the gluon propagator 
this has been done in more detail by other authors (see 

e.g. HI). 

Keeping first the lattice spacing fixed we have found 
that both the ghost and gluon dressing functions calcu- 
lated at the same physical momentum q 2 decrease as the 
lattice size L 4 is increased. This is illustrated for vari- 
ous momenta in Fig. [3J There both dressings functions 
versus the physical momentum are shown for different 
lattice sizes at either (3 = 5.8 or (3 — 6.2. Note, in this 
figure we have not dropped data with vanishing momen- 
tum components fc M to emphasize the influence of a finite 
volume on those (low) momenta. We also show data from 
simulations on a 8 4 and 12 4 lattice. One clearly sees that 
the lower the momenta the larger the effect due to a fi- 
nite volume. In comparison with (3 = 5.8 this is even 
more drastic at (3 = 6.2. At this (3 the lattice spacing is 
about a = 0.068 fm. Thus the largest volume considered 
at j3 = 6.2 is about (1.6 fm) 4 , which is even smaller than 
the physical volume of a 16 4 lattice at j3 = 5.8. 

Altogether we can state that for both dressing func- 
tions finite volume effects are clearly visible at volumes 
smaller than (2.2 fm) 4 , which corresponds to a 16 4 lattice 
at (3 = 5.8. The effect grows with decreasing momenta or 
decreasing lattice size (see the right panels in Fig. 01 . At 
larger volumes, however, the data for q > 1 GeV coincide 
within errors for the different lattice sizes (left panels). 
Even for q < 1 GeV we cannot resolve finite volume ef- 
fects for both dressing functions based on the data related 
to a 24 4 and 32 4 lattice. 

Based on our chosen values for (3 and the lattice sizes 
we can select equal physical volumes only approximately. 
Hence also the physical momenta are only approximately 
the same if the ghost and gluon dressing functions are 
compared at different (3, e.g. at different lattice spacings. 
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FIG. 3: The ghost (upper panels) and gluon (lower panels) 
dressing functions are shown for different lattice sizes as func- 
tions of momentum q 2 . The left panels show data at (3 = 5.8 
and the right ones at (3 — 6.2. Only data on be gauge copies 
are shown here. The lines are drawn to guide the eye. 
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1.0 

q 2 [GeV 2 ] 
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1.0 
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FIG. 4: The ghost (left) and the gluon (right) dressing func- 
tions at (approximately) fixed physical volume Vi ~(2.2 fm) 4 
are shown as functions of momentum q 2 . The data at (3 — 5.8 
(6.0) correspond to a lattice spacing of about a=0.136 fm 
(0.093 fm). Only data on be gauge copies are shown here. 



Therefore, it is difficult to analyze the systematic effect 
of changing a if for both dressing functions small varia- 
tions in q 2 are hidden. Consequently, in Fig. ^ we show 
the data for the ghost and gluon dressing functions ap- 
proximately at the same physical volume V » (2.2/m) 4 
for different a, albeit as functions of q 2 . This allows us 
to disentangle by eyes a change of the data due to vary- 
ing a from the natural dependence of the propagators 
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FIG. 5: The time histories of the ghost propagator calculated on first (fc) and best (be) gauge copies on a 12 4 , 16 4 and 
24 4 lattice at j3 = 6.2 and 24 4 lattice at j3 = 5.8. From left to right the corresponding runs listed in Table |H are No. 2, No.3, 
No.4/No.l0 and No.l/No.8. The upper and lower panels show data for the lowest momentum realization k — (1,0,0,0) and 
k — (0, 1, 0, 0), respectively. 



on q 2 . Inspecting Fig. 0] one concludes that the gluon 
dressing function at the same physical momentum and 
volume increases with decreasing lattice spacing. A sim- 
ilar effect (beyond error bars) we cannot report for the 
ghost dressing function. 



IV. THE PROBLEM OF EXCEPTIONAL 
CONFIGURATIONS 

We turn now to a peculiarity of the ghost propaga- 
tor at larger (3 which has also been observed by some 
of us in an earlier SU{2) study [2^|. While inspecting 
our data we found, though rarely, that there are outliers 
in the Monte Carlo time histories of the ghost propaga- 
tor at lowest momentum. Those outliers are not equally 
distributed around the average value, but are rather sig- 
nificantly larger than this. 

In Fig. we present time histories of the ghost prop- 
agator G(fc) measured of fc and be gauge copies for 
two smallest momentum realizations k ~ (1,0,0,0) and 
k = (0,1,0,0), separately. From left to right the panels 
are ordered in increasing order of the lattice sizes 12 , 
16 4 and 24 4 at j3 = 6.2 and 24 4 at (3 = 5.8. 

As can be seen from this figure in the majority ex- 
treme spikes are reduced (or even not seen) when the 
ghost propagator could be measured on a better gauge 
copy (be) for a particular configuration. Furthermore, it 
is obvious that the exceptionality of a given gauge copy is 
exhibited not simultaneously for all different realizations 
of the lowest momentum. Consequently, to reduce the 
impact of such large values on the average ghost propa- 
gator one should better average over all momentum real- 
izations giving rise to the same momentum q. This has 
been done for the results shown in Fig.^andElat least for 
the lowest momentum at (3 = 6.2. However, compared 
to the gluon propagator it takes much more computing 
time to determine the ghost propagator for all its differ- 



ent realizations of momentum q(k). 

In addition, we have tried to find a correlation of such 
outliers in the history of the ghost propagator with other 
quantities measured in our simulations. For example we 
have checked whether there is a correlation between the 
values of the ghost propagator G(k) as they appear in 
Fig. [S] with low- lying eigenvalues and eigenvectors of the 
F-P operator. They are apparent in the contribution of 
the lowest 10 F-P eigenmodes to the ghost propagator at 
this particular fc. This we shall present in a subsequent 
publication |3^| where we shall discuss the spectral prop- 
erties of the F-P operator and its relation to the ghost 
propagator. 



V. THE RUNNING COUPLING 

We shall now focus on the running coupling a s (q 2 ) as 
defined in Eq. Q where s§/(4tt) = 3/(2tt/3) for SU(3). 
Given the raw data for the gluon and ghost dressing func- 
tions on be gauge copies the average ZQ(q 2 )Zjj(q 2 ) and 
its error have been estimated using the bootstrap method 
with drawing 500 random samples. Since the ghost- 
ghost-gluon- vertex renormalization function Z\ has been 
set to one, there is an overall normalization factor which 
has been fixed by fitting the data for q 2 > q 2 to 
the well-known perturbative results of the running cou- 
pling a 2 -ioop at 2-loop order (see also 0). Defining 
x = <Z 2 /A2_ loop , the 2-loop running coupling is given by 



a2-ioop(z) = 



4tt J i 2/3 1 ln(ln: 
A) In 1 



Pi 



In a; 



(14) 



The (3- function coefficients are (3^ — 11 and (3\ = 51 
for the SU(3) gauge group and are independent of the 
renormalization prescription. The value of A 2 _i op has 
been fixed by the same fit. The lower bound q 2 has been 
chosen such that an optimal value for % 2 /dof has been 
achieved. 
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FIG. 6: The running coupling a 3 (q) as a function of momen- 
tum q 2 determined on be gauge copies. 



The results are shown in Fig. There also the 1-loop 
contribution is shown where we used the same lower 
bound q%. The best fit of the 2-loop expression to the 
data gives A 2 _ loop =0.88(7) GeV (x 2 /dof = 0.96), while 
Ai-ioo P =0.64(7) is obtained (x 2 /dof = 1.05) using just 
the 1-loop part. For q 2 we used q 2 — 30 GeV 2 . The 
value for A 2 -i 00 p is similar within errors to the SU(2) 
result given in Ref. 

Approaching the infrared limit in Fig.|S|one clearly sees 
a running coupling a s (q 2 ) increasing for q 2 > 0.4 GeV 2 . 
However, after passing a maximum at q 2 « 0.4 GeV 2 
a s (q 2 ) decreases again. Such turnover is i n ag reement 
with DSE results obtained on a torus [22l I23I |24| . There- 
fore, one can argue that this behavior is a finite lattice 
effect although we cannot resolve a difference between 
the different lattice sizes used. A similar infrared behav- 
ior for the running cou plin g has also been observed in 
different lattice studies [2(], El ■ But opposed to [2(| the 
existence of a turnover is independent on the choice of 
Gribov copy, since qualitatively, we have found the same 
behavior for a s (q 2 ) calculated on f c gauge copies. 

For completeness we mention that running couplings 
decreasing in the infrared have also been found in lattice 
studies of the 3-gluon vertex 0113 an d the quark-gluon 
vertex [36|. 

Apart from the finite volume argument given above 
to explain such a behavior, which prevents us from see- 
ing the limit a s (0) ^ mentioned in the Introduction, 
one could also put into question whether one can really 
set Zi(q) = 1 at lower momenta. A recent investigation 
dedicated to the ghost-ghost-gluon-vertex renormaliza- 
tion function Zi(q) for the case of SU(2) p| supports 
that Zi(q) w 1 at least for q > 1 GeV. 



VI. CONCLUSIONS 

We have reported on a numerical study of the gluon 
and ghost propagators in Landau gauge using several lat- 
tice sizes at (3 — 5.8, 6.0 and 6.2. Studying the depen- 
dence on the choice of Gribov copies, it turns out that for 
the gluon propagator the effect of Gribov copies stays in- 
side numerical uncertainty, while the impact on the ghost 
propagator increases as the momentum or (3 is decreased. 
However, there are indications that the influence of Gri- 
bov copies decreases as the physical volume is increased. 
This is at least expected in the light of Ref. J5| . There it is 
argued that in the continuum expectation values of corre- 
lation functions (A(xi) . . . A(x n )} over the fundamental 
modular region A are equal to those over the Gribov re- 
gion Q, since functional integrals are dominated by the 
common boundary of A and 17. Thus Gribov copies inside 
f2 should not affect expectation values in the continuum. 

While the effect of the Gribov ambiguity on the ghost 
propagator becomes smaller with increasing /3, exception- 
ally large values appear in the history of the ghost prop- 
agator in agreement with what has been observed first 
in reference |26J. These outliers we have not seen simul- 
taneously for all lattice momenta k realizing the same 
lowest momentum q(k). However, they are apparent in 
the contribution of the lowest 10 F-P eigenmodes to the 
ghost propagator at this particular k [33| . Therefore it is 
good practice to measure the ghost propagator for more 
than one k with the same momentum q(k), in order to re- 
duce the systematic errors coming from such exceptional 
values. 

We have studied the effects of finite volume on the 
one hand and of finite lattice spacing on the other. The 
first ones are found to be essential for volumes smaller 
than (2.2 fm) 4 at the same j3 whereas the discretization 
effects at the same volume are modest. Our available 
data did not allow us to extend this analysis to physi- 
cal momenta below 1 GeV, where the Gribov ambiguity 
shows up and where a similar separation of finite volume 
and discretization effects would be desirable. However, 
we could observe from Fig. [21 that enlarging the volume 
by decreasing (3 leads to a reduction of the systematic 
Gribov effect in the ghost propagator. 

The dressing functions, Zd and Zq, of the gluon and 
ghost propagators have allowed us to estimate the be- 
havior of a running coupling a s (q 2 ) in a momentum sub- 
traction scheme. Going from larger momenta q 2 to lower 
ones a s (q 2 ) is steadily increasing until q 2 « 0.4 GeV 2 . 
For q 2 < 0.3 GeV 2 , however, a s (q 2 ) is decreasing. A 
decreasing running coupling at low momenta is in qual- 
itative agreement with recent DSE results obtained on 
a torus |22l . l23l l24j . Therefore one might conclude that 
the decrease is due to finite lattice volumes we used. It 
makes it questionable whether lattice simulations in near 
future can confirm the predicted infrared behavior of the 
gluon and ghost dressing functions with related expo- 
nents kd = 2kg- 
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APPENDIX A: SPEEDING UP THE INVERSION 
OF THE F-P OPERATOR 

For the solution of the linear system M<f> = ip c with 
symmetric matrix M, the conjugate gradient (CG) al- 
gorithm is the method of choice. Its convergence rate 
depends on the condition number, the ratio of largest to 
lowest eigenvalue of M. When all U Xjll — t obviously the 
F-P operator is minus the Laplacian A with a diagonal 
color substructure. Thus instead of solving M <j) = ip c one 
rather solves the transformed system 

[MA" 1 ] (A0) = tf c 

In this way the condition number is reduced, however, the 
price to pay is one extra matrix multiplication by A -1 
per iteration cycle. In terms of CPU time this should be 
more than compensated by the reduction of iterations. 

The pre-conditioned CG algorithm (PCG) can be de- 
scribed as follows: 

initialize: 

7 (0) = ( p -(o) jr -(o)) 
start do loop: k = 0, 1, . . . 

g(k+l) = A -l-(fe+l) 
7 (fc+l) = (^(fc+l) 5 f(fe+l)) 

if (^ k+1%) < e) exit do loop 

end do loop 

Here (•, •) denotes the scalar product. 

To perform the additional matrix multiplication with 
A -1 we used two fast Fourier transformations J 7 , due to 
(-A)" 1 = T~ Y q~ 2 (k) T. The performance we achieved 
is presented in Table [Q] We conclude that on larger 
lattice sizes the reduction of iterations is about 70-75%, 
while the resulting reduction of CPU time depends on 



the lattice size. This is because we are using the fast 
Fourier transformations in a parallel CPU environment. 
If the ratio of used processors to the lattice size is small 
(see e.g. the data for 32 4 lattice at this table), almost 
the same reductions of CPU time as for the number of 
iterations is achieved. 

Further improvement may be achieved by using 
the multigrid Poisson solver to solve Az^ = f^ k ' . 
This method is supposed to perform better on 
parallel machines. Perhaps a further improve- 
ment is possible by using as pre-conditioning matrix 
M -1 = — A -1 — A _1 AfiA _1 + . . . which is an approxi- 
mation of the F-P operator M = — A + Mi to a given 
order (see also [2j). However, the larger the order, 
the more matrix multiplications per iteration cycle are 
required. This may reduce the overall performance. We 
have not checked so far which is the optimal order. 





CG 




PCG 




speed up 


lattice 


iter 


CPU [sec] 


iter 


CPU [sec] 


iter CPU [sec] 


8 4 


1400 


3.7 


570 


2.4 


60% 35% 


16 4 


3900 


240 


1050 


130 


73% 46% 


32 4 


9900 


13400 


2250 


3900 


77% 71% 



TABLE II: The average number of iterations and CPU time 
per processor (PE) using the CG and PCG algorithm to in- 
vert the F-P operator are given. All inversions have been 
performed at f3 — 5.8 with source S bc exp(27ri k ■ y) where 
k = (1,0,0,0). To compare the different lattice sizes 4 PEs 
have always been used. 
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